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We analyze canonical determinants, i.e., grand canonical determinants projected to a fixed net 
quark number. The canonical determinants are the coefficients in a fugacity expansion of the 
grand canonical determinant and we evaluate them as the Fourier moments of the grand canonical 
determinant with respect to imaginary chemical potential, using a dimensional reduction technique. 
The analysis is done for two mass-degenerate flavors of Wilson fermions at several temperatures 
below and above the confinement/deconfinement crossover. We discuss various properties of the 
canonical determinants and analyse the convergence of the fugacity series for different temperatures. 



PACS numbers: ll.15.Ha 

I. INTRODUCTION 

QCD with finite density is an important topic of both, 
experimental and theoretical studies. In particular one 
would like to understand the various transitions to dif- 
ferent states of matter that are conjectured for the QCD 
phase diagram. Exploring these phase transitions is 
clearly a non-perturbative problem and non-perturbative 
methods must be applied for such studies. In princi- 
ple the lattice QCD formulation provides a suitable non- 
perturbative framework, and as long as the chemical po- 
tential is zero this approach may be used to obtain reli- 
able quantitative results. 

However, if one works at finite density a severe techni- 
cal problem emerges: The fermion determinant is com- 
plex at non-zero chemical potential n and cannot be used 
as a probability weight in a Monte Carlo calculation. Var- 
ious reweighting and expansion techniques have been ex- 
plored in recent years and the reviews at the annual Lat- 
tice conference provide a glimpse at the corresponding 
developments [Ij . An interesting option is to work not in 
the grand canonical approach with a chemical potential 
jjL, but instead to use canonical partition sums with fixed 
net quark number q. Several studies of this alternative 
perspective can be found in the literature [2] - [in]- So 
far the focus was mainly either on reweighting techniques 
or full canonical simulations at a fixed quark number. 

Here we explore the perspectives of a fugacity expan- 
sion of the grand canonical ensemble. The expansion 
coefficients are fermion determinants projected to a fixed 
net quark number, so-called canonical determinants. The 
expansion is a Laurent series in the fugacity parameter 
e^^ {(3 is the inverse temperature) which has properties 
different from the more conventional Taylor expansion in 
fif3. However, the possible application to lattice QCD at 
finite density is the same for both expansions: One com- 
putes the expansion coefficients in simulations at = 



and uses the series to explore QCD at /i > 0. 

In this article we discuss the evaluation and the prop- 
erties of canonical determinants. They are computed as 
Fourier moments of the fermion determinant at imagi- 
nary chemical potential fi/3 — itp, which we evaluate effi- 
ciently using a recently proposed dimensional reduction 
formula [8]. We study the distribution properties of the 
canonical determinants and analyze the convergence of 
the fugacity expansion for various temperatures in the 
confined and deconfined phases. This analysis is done 
for a fixed numerical effort, i.e., a fixed number of 256 
sampling points of cp in the interval (— tt, tt] used for the 
evaluation of the Fourier transformation to obtain the 
canonical determinants. 



II. CANONICAL DETERMINANTS AND 
THEIR EVALUATION 

The starting point for our study of canonical deter- 
minants and the fugacity expansion is Wilson's lattice 
Dirac operator with chemical potential fj, on lat- 

tices of size X Nt- The lattice spacing will be denoted 
by a, such that /3 = aNt- For two mass-degenerate flavors 
the corresponding grand canonical partition function is 
given by 

Z{^,) = J !?[[/] e-^^[^l det[D{^i)]\ (1) 

where Sa[U] is the action for the gauge fields U, 2?[?7] the 
path integral measure and det[£'(/i)] the fermion deter- 
minant for chemical potential /i, which in the following 
will be referred to as grand canonical determinant. The 
grand canonical determinant can be expanded in a finite 
fugacity series, 

det[D{f,)] = Y.^^^'D^'^ ' (2) 

9 
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where the sum runs over integer valued quark numbers 
q € [—6N^,+6Ng]. The expansion coefficients D^^' are 
tlic canonical determinants and may be obtained us- 
ing Fourier transformation with respect to an imaginary 
chemical potential, 

^ ^ f dip e-"?'^ det[D(^/3 = i(p)] . (3) 

While for real chemical potential /i the grand canonical 
determinant is complex, for imaginary chemical poten- 
tial the determinant det[_D(/i/3 = iip)] is real. One may 
explore the generalized 75-hermiticity relation for Wilson 
fermions, "f^D{ii)"f^ — D{~y)\ to estabhsh the relation 
= (£)(?))* between canonical determinants with 
positive and negative quark numbers q. For vanishing 
quark number q = Q the canonical determinant is real, 
i.e., i:>(0) G R. 

In principle Eq. ([s]) is a very elegant expression for the 
canonical determinants D'^'^\ However, the Fourier inte- 
gral has to be evaluated numerically and the evaluation 
of the grand canonical determinant det[Z3(/i/3 — iLp)] for 
the necessary (/^-values is computationally rather costly. 
To alleviate the problem we use a domain decomposition 
technique |5] where the grand canonical fermion determi- 
nant may be rewritten exactly in a dimensionally reduced 
form, 

det[i:i(^)] = Aodet [I-Hq- e^^H - e^^'^ff ^] , (4) 

which holds both for real and imaginary chemical poten- 
tial. Here Aq is a real factor which depends only on the 
background gauge field configuration but is independent 
of the chemical potential /i. Hq = Hq, and H are ma- 
trices that are built from propagators on sub-domains of 
the lattice and live on only a single time slice (see [8] 
for details). Thus the determinant in Q is dimension- 
ally reduced, i.e., the determinant is taken over a matrix 
with 12A^^ rows and columns, where the factor 12 comes 
from the color and Dirac indices. The terms Hq and H 
can be completely stored in memory and are then used 
many times for the evaluation of det[_D(/i/3 = i(f)] at the 
necessary values of ip. Due to the reduced size of the ma- 
trix where the determinant has to be evaluated, we gain 
a factor of 0{N^), which here, at Nt = 4, corresponds to 
a speedup by a factor of 64. The Fourier integral (|3| is 
evaluated numerically with standard techniques ;11J us- 
ing 256 values of ip in the interval (— 7r,7r]. We remark 
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TABLE I: Parameters of our numerical simulation. 



at this point that also alternative dimensional reduction 
formulas with different properties were proposed [21 1121 - 

[g. 

Our numerical tests are done on dynamical configu- 
rations with two flavors of Wilson fermions at = 
generated with the MILC code [TS]. In this exploratory 
study of the canonical determinants and the fugacity ex- 
pansion we are here limited to a single lattice volume, 
8"^ X 4. The parameters (inverse gauge coupling Pgauge 
and hopping parameter k) and the values for the cor- 
responding lattice spacing a and the temperature T are 
taken from [5] and listed in Table |lj Due to the small 
lattice the pion mass has to be kept rather large - it is 
close to 950 MeV for all our ensembles. For each set of 
parameters we use a statistics of 200 configurations. All 
errors we show are statistical errors determined with the 
jackknife method. In addition to the parameters of the 
simulation, in Table [T] we also list the cutoff parameter 
qcut used in the fugacity expansion which we discuss in 
more detail in Section IV. 

In the range of temperatures we consider, the lowest 
three values of T are in the confined phase, the crossover 
into the deconfined phase is at T ^ 165 MeV for our 
simulation parameters, and the two largest temperatures 
are in the deconfined phase. 

As a first result we now have a look at the average 
(logj^g det[_D(/i/3 — iip)]) of the logarithm of the grand 
canonical determinant as a function of the imaginary 
chemical potential parameter p. This is interesting be- 
cause the Fourier moments with respect to Lp are the 
canonical determinants we want to study and the Lp- 
dependence of det[£'(/i/3 — ip)] already provides a first 
insight into their properties. (..) is the average over 
gauge configurations generated for two mass degenerate 
flavors of Wilson fermions at /i = as detailed above. 
We study the (/s-dependence of (log^Q det[D(/i/? — iip)]) 
using all 256 values of p in the interval (— 7r,7r]. The 
integrands det[_D(/i/3 — ip)] are computed with the di- 
mensional reduction formula Q as outlined. We show 
the corresponding results for all our temperature values 
in Fig.[l} 

The figure clearly demonstrates that in the confined 
phase the grand canonical determinant shows only a 
weak dependence on ip^ while above the crossover a 
pronounced variation with ip develops. This behavior 
has a simple and well known interpretation: An imag- 
inary chemical potential — iip is equivalent to an 
additional temporal boundary condition for the fermion 
fields, '(p{x, Nt + 1) = —e^'^4'{x, 1), where the minus sign is 
the usual anti-periodic temporal boundary condition for 
fermions. While in the confined phase the correlations 
are short ranged and do not give rise to a dependence 
on the boundary angle ip, in the deconfined phase cor- 
relations over distances larger than the temporal extent 
of the lattice generate the non-trivial response to chang- 
ing ip. The canonical determinants Z)''^ are the Fourier 
moments of the (^-dependence and the first analysis in 
Fig. [1] already demonstrates that in the deconfined phase 
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FIG. 1: Average of the logarithm of the grand canonical 
fermion determinant as a function of the boundary phase ip. 
We compare the results for different temperatures from 144 
MeV (bottom curve) to 211 MeV (top). 



the higher modes, i.e., Z)'^''^ with larger values of \q\ will 
play a more prominent role. 



III. PROPERTIES OF THE CANONICAL 
DETERMINANTS 

Having discussed the canonical determinants D^'^ and 
the setting of our calculation we now come to analyzing 
properties of the In Fig. [2] we show scatter plots 

in the complex plane for D^'^^ / det[D{ii ~ 0)], i.e., the 
canonical determinants normalized by the grand canoni- 
cal determinant at /i = 0. This normalization is such that 
when one sums over q all data points / det[D(/x = 0)] 
in the complex plane for a single configuration the result 
is 1. In the plot we show results for g = 0, 1, 2, 3 and com- 
pare T = 144 MeV (top row of plots in Fig. [2]), T = 164 
MeV (middle) and T = 211 MeV (bottom). Note that 
different axis scales are used. 

As already discussed, for g = the canonical determi- 
nants d(9=o) must be real, a fact that is obvious in the 
figures. The corresponding values scatter on the real axis 
for all temperatures. For g > the points spread out 
also in the imaginary direction. In the confined phase 
(T = 144 MeV) the points for the higher Fourier modes 
L*'^' then quickly approach the origin of the complex 
plane when increasing q. The situation is different in the 
deconfined phase {T = 211 MeV), where we observe that 
for increasing q the values of Z?^''^ do not move toward 
the origin as quickly as for T = 144 MeV. For T = 164 
MeV (middle row of plots in Fig. [2]) one finds an in- 
termediate behavior. The relative size of the canonical 
determinants I?^'^-' at the different temperatures refiects 
the observation we have already made in the discussion of 



Fig. [T] For temperatures above the crossover transition 
higher Fourier components D^'^^ are necessary to resolve 
the strong dependence of the grand canonical determi- 
nant det[D(/i/3 = lip)] on the boundary phase tp. 

We conclude the discussion of Fig. [2] with noting, that 
the scatter plots of D'^'^-' for the smaller values of q show a 
curious pattern, in particular the q = I data for T — 144 
and 164 MeV: They arrange in two oblong structures that 
seem to have a preferred angle relative to the real axis. 
A similar pattern at temperatures near the crossover was 
observed also in canonical determinants for staggered 
fermions |16| . At the moment we do not understand the 
origin of this pattern, but it might be related to the prop- 
erties of the D^''^ under center transformations which 
rotate them according to D^'^'> ^'^ where 
z e {1, e*^'^/'^, e"'^'^/'^} is the center element multiphed 
to all temporal gauge links in a fixed time slice. For pure 
gauge theory the spontaneous breaking of center symme- 
try is intimately related to the deconfinement transition, 
and it could be that the pattern we observe is a remnant 
of this mechanism. This speculation is further supported 
by the observation, that the phase of the Polyakov loop 
(the order parameter for confinement in pure gauge the- 
ory) is strongly correlated with the phases of the D'-'^^ 
|16) . Also in an analysis of the canonical determinants 
of the quenched case [HI [TO] it was demonstrated that 
the D'-'^^ in the broken phase very cleanly map the center 
orientation of the underlying center sector. 

The change of the relative weight of the Z?^'^-' with tem- 
perature is manifest also in Fig. [3j where we study the 
size distribution (|_D'^' l/D^")) of the canonical determi- 
nants normalized to the trivial q = sector as a function 
of q and again compare results for T = 144 MeV (top), 
T = 164 MeV (middle) and T = 211 MeV (bottom). 
The plots show that the distribution roughly follows a 
Gaussian centered around the dominant q = sector. 
While in the confined phase the distribution is rather 
narrow, in the deconfined phase it widens and reflects the 
fact that at higher temperatures also sectors with larger 
quark numbers become accessible. When one turns on 
a chemical potential /i > the canonical determinants 
are multiplied with powers of the fugacity factor, i.e., 
e^'5«D(«), which has the effect of shifting the 
center of the distribution towards larger quark numbers. 

Finally, in Fig.jljwe have a look at the phase 6*'^-' of the 
canonical determinants D^'^\ More precisely we look at 
(e^2e<''^ = (dW/dC-?)) (remember that = D^?)*). 

In Fig. [4] we show (D^'^^/f^"')) as a function of q for 
different values of the temperature. 

Below the crossover temperature the results for 
(D(«V^^"'^) drop to rather quickly with mcreasmg g, 
i.e., the phases of the Z?^'^-' for q above roughly q — 5 fluc- 
tuate strongly and average to a very small number. For 
the largest temperatures the decrease of {D^'^^ / D^^''^) 
with q is slower and the phases have a sizable expectation 
value up to g ~ 25 (we stress that these are statements 
specific to the parameters of our calculation, in particular 
the volume used). 
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FIG. 2; Scatter plot of the canonical determinants Z)''-*/ det[_D(/i = 0)] in the complex plane at T = 144 MeV (top row of 
plots), T = 164 MeV (middle) and T = 211 MeV (bottom) for q = Q,l, 2, 3. The canonical determinants are normalized by the 
grand canonical determinant at /i = 0. Note that we use different scales for the different temperatures. 



The decrease of the (D^'V-D^"''^) with increasing q 
also sheds light on how the fermion sign problem may 
be viewed in the fugacity expansion: When increasing 
^ the powers e^*^' of the fugacity factors put a larger 
weight on the canonical determinants ZJ^"?) with higher 
quark numbers, which in turn arc the ones that are sup- 
pressed due to the fluctuations of their phase. The fact 
that for higher temperatures the sign problem is milder 
is clearly visible from the slower decrease of (Z)'') / D^~'^^) 
for larger temperatures. 



IV. CONVERGENCE PROPERTIES OF THE 
FUGACITY SERIES 

After the study of the canonical determinants _D^^', 
we now can start to analyze the fugacity series ([2]). In a 
numerical study it is necessary to truncate the series and 



we denote the truncated fugacity expansion by 



E 



(5) 



where qcut denotes the highest (anti-) quark sector we 
take into account in the truncated series. 

The central question we want to study in this section is 
how well the truncated series S{ii)q^^^^ approximates the 
full grand canonical determinant det[£'(/i)] and how the 
quality of the approximation depends on the parameters 
fi and T. This is a question which we address both for 
individual configurations and for the gauge average. We 
stress at this point, that for a finite volume and non-zero 
temperature, the representation of det[£'(^)] by S{^)q^^^ 
would be exact with qcut = 6A^^ if all coefficients D^'^^ 
were known with arbitrary precision. Thus as long as the 
lattice is finite such that no singularities can emerge, an 
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FIG. 3: Distribution of {\D'''^\/D^°^) as a function of the 
quark number q. We compare three temperatures: T = 144 
MeV (top plot), T = 164 MeV (middle) and T = 211 MeV 
(bottom). 



insufficient representation of the grand canonical deter- 
minant is only due to limited accuracy in the numerical 
evaluation of the Fourier integrals ^ . 

We assess the quahty of the approximation by consid- 
ering the relative error 



'S'(M)<7c„t 



(6) 



For the study of A we also evaluated the grand canonical 
determinant det[Z)(/x)] for a few values of real chemical 
potential, again applying the dimensional reduction for- 
mula These values are then used for the assessment 
of the fugacity expansion using A. 

Inspecting the fugacity series (|2]) or ([5| one observes 
that for increasing q there is competition of two terms: 
The powers of the fugacity, i.e., the factors e'^^', are 
terms that grow exponentially with q and their rate of 




FIG. 4: Expectation value of the phase of the canonical de- 
terminants, (e"***'') = (D^'V-D'""') as a function of q. We 
compare the results for different temperatures. 



growth is determined by the product fi/3. The expo- 
nential growth with q is compensated by the quick de- 
crease of the canonical determinants D^'^\ This decrease 
is roughly Gaussian as we demonstrated in Fig. [3] and 
thus dominates the increasing function e''*^*. However, 
this mechanism can be spoiled by numerical instabilities: 
If the higher coefficients D*^'-' are not known with perfect 
accuracy (as will always be the case when they are eval- 
uated via numerical Fourier transformation), their fluc- 
tuations will be amplified by the factors e'^^'. This ef- 
fect is manifest in the behavior of S{fj,)q^^^ when q^ut is 
varied: First S{fi)q^^^ quickly saturates as a function of 
qcut, but for too large values, e.g., already qcut ^ 30 for 
some of our parameter values, the series S{fi)q^^^ starts 
to diverge. Fortunately the D*^*^ decrease very rapidly: 
Fig. |3] shows that also for our highest temperature en- 
semble all D'^'^^ with l^l > 20 essentially vanish, and the 
series S{^)q^^^ can be truncated at small qcut- The op- 
timal values for q^ut were determined using the relative 
error A and we list them in Table HI 

In Fig. [5] we show how the relative error A of the fu- 
gacity expansion depends on the parameters /i/3 and T. 
For all 200 configurations in each of our ensembles we 
compute the canonical determinants as well as the 

canonical determinant det[D{^)] for several real values of 
fi. For these values of we can evaluate A. 

Fig. [5] shows A as a function of the dimensionless com- 
bination /z/3, and we compare the results for all available 
temperatures. For each ensemble we show A for the con- 
figuration where A is largest ("worst"), for the config- 
uration where A is smallest ("best"), for the median of 
the 200 configurations, as well as the average of A over 
all configurations. 

From Fig. [5] it is obvious that for the lower values of T 
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FIG. 5; Relative error A (see the text for its definition) of the fugacity expansion as a function of at different temperatures. 
We show the error of the best configuration, the worst configuration, the average over all configuration as well as the median. 



the fugacity expansion has better convergence properties: 
At r = 144 MeV all four categories (best, worst, median 
and average) have a relative error A smaller than 1 % 
for values of the chemical potential up to /i/3 — 0.8 and 
the average error over all configurations remains below 
1% even up to = 1.1. For T ~ 211 MeV an error 
smaller than 1 % can be maintained only up to = 0.4. 
The reason is that here the distribution of the Z?''^ is 
much wider (compare Fig. [s]), and the D^'^^ for larger 
values of Iq'I contribute significantly in the fugacity sum 
already at not too large /Lt/3. Since these higher Fourier 
modes are only known with less relative accuracy, the 
fugacity expansion at the highest temperatures breaks 
down already at smaller chemical potentials. 

We stress that our convergence analysis is specific for 
the setting we use here, i.e., lattice size 8'^ x 4, parameters 
as listed in Table |l] and 256 values of in the evaluation 
of the Fourier integrals ([3|. In particular increasing the 
number of sampling points for in the Fourier integral 
will allow for a higher precision of the D*^''-* which leads to 
better convergence properties of the fugacity sum, such 
that higher values of [iji can be reached. 

An interesting question is how the number of D^'^-' 
that need to be evaluated depends on the volume of the 
box. This number is roughly proportional to the width 
of the distribution of the D'-'-' as displayed in Fig. [sj In 
the next section we will show that this width is related 



to the quark number susceptibility Xq- This is an ex- 
tensive quantity, i.e., it grows with the spatial volume 
V — (aNsf' . We thus conclude that the number of D'^'^'^ 
that contribute significantly to the fugacity sum is pro- 
portional to the spatial volume Y . In turn this means 
that also gcut needs to be increased linearly in Y . 



V. QUARK NUMBER DENSITY AND QUARK 
NUMBER SUSCEPTIBILITY 

Let us finally discuss an exploratory calculation of the 
quark number density and the quark number susceptibil- 
ity based on canonical determinants and the fugacity ex- 
pansion. Combining the general expression for the grand 
canonical partition sum Z(/i) from Eq. ([T]) with the fu- 
gacity expansion ^ we find 

Z(/i) - y2?[;7]e-^«[^l^^e'^^i?(')^ (7) 

where in the second step we inserted 1 = 
det[£'(0)]^/ det[£)(0)]^ to write the whole expression as 
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a vacuum expectation value (up to normalization) in the 
grand canonical ensemble at /i = where we perform 
our simulation. The definitions of the quark number 
density Uq and the quark number susceptibility Xq are 



r3 



2^2 



/33 lnZ(/i) 



(8) 



Both observables are intensive quantities after the nor- 
malization with the 3- volume V and are made dimension- 
less using suitable powers of T. For the necessary first 
and second derivatives one finds 
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where we introduced the moments M„ , n 
fugacity series as 
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(10) 



Finally we express the normalization factor 1/Z{^) with 
the moment Mo and obtain (use jS^/V = (Nt/Ng)^) 



MqMi ) 
(Mo)2) ' 
{MoAh + jAhf] 
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((Mo)^; 



(11) 
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In the final result ( 11 ) both observables are expressed in 
terms of vacuum expectation values of moments of the 
fugacity series. These vacuum expectation values (..) are 
computed with two flavors of Wilson fermions at fi — 0. 

In Fig. 6 we show our results (symbols) for Uq/T^ as 
function of the dimensionless combination /i/3 for three 
different temperatures. We compare the results to the 
outcome of the Taylor expansion presented in [T71 [TH] . 
More precisely we used the terms up to fourth order from 
[T5] and interpolated the Taylor coefficients to the tem- 
peratures used in our study. However, this comparison 
should be viewed with caution: The calculation in (18j is 
done on considerably larger lattices, has a statistics that 
is by a factor 20 larger than the statistics available in 
our exploratory study, and also is based on the staggered 
formulation. Nevertheless we find reasonable agreement 
for our results at these temperatures and range of 

We also attempt a comparison of our results for the 
quark number susceptibility Xg at /i = to the results 
from the Taylor expansion. The corresponding numbers 
are listed in Table [TT| (again interpolated to our values of 
T). We find that up to temperatures near the crossover 
the agreement of the two approaches is reasonable. Given 
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FIG. 6: Quark number density as a function of jj,/3. We com- 
pare our results at three different temperatures (symbols) to 
the results from Taylor expansion (curves). 



the fact that the calculation [18] and our exploratory 
study differ considerable in statistics and volume, this 
agreement is satisfactory. Above the crossover we see, 
however, a rather strong discrepancy. We suspect two 
main reasons for this discrepancy: 1) The rather small 
volumes that are available for our calculation, and 2) the 
fact that for the larger values of T higher Fourier modes 
contribute significantly (the factor in the second mo- 
ment M2 of Eq. ( 10 ) enhances them further) which in 
the current setting of 256 sampling points for the numer- 
ical evaluation of the Fourier integral are not sufficiently 
accurate. 



VI. SUMMARY AND DISCUSSION 

In this exploratory study we analyzed the canonical de- 
terminants of 2-flavor lattice QCD with Wilson fermions. 
The canonical determinants are the coefficients in the fu- 
gacity expansion of the fermion determinant and may be 
calculated as the Fourier modes with respect to imag- 
inary chemical potential. We speed up this evaluation 
by using a dimensional reduction formula for the grand 
canonical determinant. 



T 



Fugacity expansion Taylor expansion 18' 



144 MeV 
153 MeV 
157 MeV 
164 MeV 
173 MeV 
189 MeV 
211 MeV 



0.25(8) 
0.14(9) 
0.41(8) 
0.52(7) 
0.90(8) 
3.41(8) 
5.85(4) 



0.1360(44) 
0.2312(52) 
0.3190(58) 
0.5340(62) 
0.9170(66) 
1.3514(36) 
1.5240(30) 



TABLE II: Comparison of results for Xq/T^ at = 0. 
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We illustrate that a sizable dependence of the fermion 
determinant on the imaginary chemical potential sets in 
only at temperatures near the crossover temperature, 
which already shows that at the crossover an enhance- 
ment of the higher quark numbers can be expected. 
Studying the size distribution of the canonical determi- 
nants we find that their distribution is roughly Gaus- 
sian in the net quark number q with a width that starts 
to increase at the crossover. The average phases of the 
canonical determinants drop quickly with q, but the drop 
is slowed above the crossover temperature. This analysis 
sheds light on the complex phase problem from the point 
of view of the fugacity series. 

We continue with a systematical analysis of the conver- 
gence properties of the fugacity expansion by comparing 
the truncated fugacity series to an exact evaluation of 
the grand canonical determinant at several values of the 
chemical potential. It is shown that for lower tempera- 
tures we can obtain very good accuracy up to ^ 0.8, 
while above the crossover we only reach /i/3 ^ 0.4. We 
stress that these results are specific for the numerical 
effort we invest here, in particular 256 sampling points 
for the evaluation of the Fourier moments. On a finite 
lattice the fugacity expansion is a finite Laurent series, 
and in principle it is possible to compute all coefficients 
such that the series representation of the grand canonical 
determinant becomes exact. 

Finally we use the fugacity series to explore the evalua- 
tion of the quark number density and the quark number 
susceptibility through first and second moments of the 
truncated fugacity series. This is an interesting alterna- 
tive to a standard calculation where the quark number 
g is a binomial in the quark fields, q oc J d'^xip{x)j4ip{x) 
which after integrating out the fermions is related to the 
trace of the quark propagator, g c>c Tr j4D~^. In this 
form q is difficult to evaluate numerically using, e.g., 
stochastic estimators. Once the canonical determinants 
Z?^^^ are available the expressions in terms of the mo- 
ments of the fugacity expansion are considerably simpler 
and we demonstrate that the approach is compatible with 
the conventional method up to temperatures where we 
have sufficiently control over the fugacity series. 

Conceptually the fugacity expansion falls in the same 
category as the Taylor expansion: A series is used to ex- 
trapolate the information from a Monte Carlo simulation 



at /i = to finite values of the chemical potential. Thus 
both, Taylor- and fugacity expansion face the same over- 
lap problem, i.e., the /i = configurations used for eval- 
uating the expansion coefficients may have rather little 
overlap with the configurations that dominate physics at 
finite /i. This limits both expansion approaches to small 
values of fi. 

Concerning other properties, Taylor- and fugacity ex- 
pansion are different: At finite volume the fugacity series 
is a finite Laurent series in the fugacity parameter e'^'^, 
while the Taylor series is an infinite series even at finite 
volume. Thus the two series will have different proper- 
ties and one should explore which of the two approaches 
provides a better expansion around /i = 0. To study this 
question we currently compare the two series in an ef- 
fective theory for QCD, where a fiux representation free 
of the complex phase problem allows one to obtain high 
precision reference data from a Monte Carlo simulation 
19 . These will be used to study the quality of the ap- 
proximation from Taylor- and fugacity series. 

We conclude with stressing that this exploratory work 
should be considered only as a first step towards a full 
development of fugacity expansion as a reliable tool in fi- 
nite density lattice QCD. Certainly further improvement 
of the numerical techniques and more tests will be nec- 
essary to assess the full potential of the approach. 
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